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Bilayer graphene hosts valley-chiral one dimensional modes at domain walls between regions of 
different interlayer potential or stacking order. When such a channel is brought into proximity to a 
superconductor, the two electrons of a Cooper pair which tunnel into it move in opposite directions 
because they belong to different valleys related by the time-reversal symmetry. This is a kinetic 
variant of Cooper pair splitting, which requires neither Coulomb repulsion nor energy filtering but 
is enforced by the robustness of the valley isospin in the absence of atomic-scale defects. We derive 
an effective model for the guided modes in proximity to an s-wave superconductor, calculate the 
conductance carried by split and spin-entangled electron pairs, and interpret it as a result of local 
Andreev reflection processes, whereas crossed Andreev reflection is absent. 

PACS numbers: 72.80.Vp,74.45.+c,03.65.Ud 


Creating mobile nonlocal spin-entangled electrons in a 
transport experiment with the help of superconductor- 
normal junctions has attracted a lot of attention in theory 
m and experiment [5HI1] because the spin degree of 
freedom of the electron could serve as a solid-state qubit 
m In the existing experiments, the envisaged process 
where a Cooper pair is split over two normal leads is 
crossed Andreev reflection (CAR) [161 IE] j which is en¬ 
hanced by the repulsive electron-electron interaction on 
two quantum dots weakly coupled to the superconduc¬ 
tor [T] or by energy filtering H HHj. The basic mech¬ 
anism of these entanglers is not very sensitive to the 
specific material used, i.e., the underlying band struc¬ 
ture. It has been shown that characteristic features of 
new materials exhibiting Dirac-cones like graphene or 
topological insulators can be useful for splitting Cooper 
pairs [TllIMl]. In these proposals, the efficiency of the 
splitting process, in the absence of interactions, relies on 
non-protected resonance conditions or the split Cooper 
pair is not spin-entangled due to spin-helicity or spin- 
polarization of the leads. Helical edge states of the quan¬ 
tum spin Hall regime have, however, been proposed to 
detect spin entanglement (3113] . 

Here, we propose to exploit the valley degree of free¬ 
dom in bilayer graphene (BG), where valley-chiral, spin- 
degenerate one-dimensional (ID) channels are formed at 
domain walls. Such domain walls can be engineered by 
switching the sign of an interlayer voltage or by revers¬ 
ing the stacking order PH US] . If brought into proximity 
to a superconductor, the pairs emitted into the chan¬ 
nel are split, i.e., two electrons propagate in different 
directions but remain spin-entangled since, as required 
by time-reversal symmetry, the two electrons forming 
the Cooper pair in the superconductor are from differ¬ 
ent valleys [3S]. As long as the valley degree of free¬ 
dom is robust, the splitting efficiency is unity, indepen¬ 
dent of resonance conditions. The device extends the 
upcoming “valleytronics” in graphene |27j to nonlocal 



FIG. 1. (a) At a domain wall (red) in BG between different 
interlayer voltages ±Vb or different stacking order a topologi¬ 
cal valley-chiral channel forms. Gooper pairs tunneling into it 
from a nearby s-wave superconductor (S) are split because the 
two electrons belong to opposite valleys K± and thus have op¬ 
posite velocities. They remain spin entangled and propagate 
to separate normal leads (N). (b) In each valley two subgap 
modes along the domain wall emerge. Energy and momentum 
conservation along the NS interface single out four points in 
the subgap spectrum at which Gooper pairs are injected. 


Einstein-Podolsky-Rosen pairs. A ID channel defined by 
opposite stacking order has recently been created exper¬ 
imentally in BG |28j with mean free paths over several 
100 nm, demonstrating weak intervalley scattering. In 
this scenario, the normal reflection of an incoming hole 
(or electron) and CAR are absent. In the limit of a weak 
proximity effect, where the normal transmission through 
the proximity region of the channel is almost perfect, 
the spin-entangled pair emission with electrons moving 
to opposite normal terminals (Fig. is equivalent to lo¬ 
cal Andreev reflection (FAR) processes, opposite to the 
normal reflection-dominated case, where CAR produces 
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entangled pairs. 

We analyze the setup of Fig. [^in two steps: first, we 
investigate the influence of the superconductor on the 
ID channel by solving a Bogoliubov-de Gennes (BdG) 
equation, and derive an effective ID model to describe 
the proximity effect in the channel. Second, we calculate 
the subgap conductance when applying a bias voltage 
between the superconductor and the channel using a rate 
equation approach. We interpret the subgap transport in 
a scattering matrix picture and show that the outgoing 
scattered state is a two-particle spin-entangled state on 
top of a filled normal-state Fermi sea with a chemical 
potential lowered by the bias voltage. To leading order 
its weight is given by the LAR amplitude. 

Model. —We consider a BG sheet with Bernal AB 
stacking in the presence of an interlayer voltage V{r) 
P5] . We model the superconductor region as BG in 
which the bands are shifted by a scalar potential U{r) 
due to doping and which has an induced s-wave pair¬ 
ing amplitude A(r). We employ the low-energy approx¬ 
imation for BG, valid at energies and (inter)layer volt¬ 
ages smaller than the interlayer hopping 71 ~ 0.3 eV. 
Without the superconductor, the valley index = ±1 
distinguishing the two iF-points K± = ±( 47 r/ 3 a, 0 ) and 
the electron spin s = ±1 = 1 , 4 - are good quantum num¬ 
bers and we write the Bogoliubov-de Gennes equation as 

^BdG = i^Xvdxdyay + - d^)a^) 

+U{r) + V{r)a^]T^ +A{r)T^, (1) 

where a = Vp/ji. The Pauli matrices Ci act in the 
pseudospin {Ai,B 2 ) space and Ti in electron-hole 
space and we set the Fermi energy Ep = 0. The 
4-component spinor is 

where we have introduced the electron = 

'^S 2 ,x«,«(’’)) components [u —>■ v) 

on the two sublattices. Excitations with energy e are 
then expanded as 7 x„,s(£) = / ' 'I'x« 

with the vector of field-operators = 

In the absence of the superconductor (A(r) = U{r) = 
0 ) and assuming the modes to propagate along the y- 
direction along a domain wall at x = 0, i.e., V{r) = 
—Vo sgn(x) with Vb > 0 , the topologically confined modes 
can be found analytically [SJ. The electron and hole 
sectors in Eq. Q decouple. In the electron sector, the 
solutions in each half space have the form = 

s(r),0)'^ where ^{r) = ^qHp^^+Pvv) with 

For any fixed energy e and momentum py there are four 
allowed values Px = Ei^/Vo — jet — Py which be¬ 
come complex when |e| < Vb, i.e. there are no propagat¬ 
ing modes in the bulk at energies below Vb. Matching the 


wavefunctions decaying away from the domain wall and 
their derivatives, one obtains the two electronic subgap 
solutions n = 1,2 in each valley, $®’"gp^(x), with the 
dispersion relation 


ai/2 I ^^0 - ^pI 

^Xv,S,Py 2 


y2V2aVb -k a2p2 ^ (3) 


with velocities opposite in the two valleys [21]. The so¬ 
lutions for the hole sector = (0, v^^^^(r))'^ 

where ^(r) at energy are obtained from 

Eqs. ([^ and Q by setting e —>• —e. 

The relevant momenta Py are close to the K points: 
taking £ ~ 0, we obtain from Eq. ([^ the momentum 
scale Py ^ \JVola, on which the K-points are located 
at ^Ix/Vdla ~ 102 for E ~ A ~ meV. The guided 
modes decay into the bulk on a length scale of x/WaJVo, 
which then is on the order of several 10 nm. This sets 
the scale of the separation between the guided mode and 
a superconductor required to obtain a proximity effect. 

Perturbation theory for superconducting pairing .— 
Assuming a superconductor/normal interface with trans¬ 
lational invariance along the y-direction, there are three 
distinct areas: in the superconductor area, x < —d, the 
pairing amplitude A(r) = A is finite and U{r) = —Us is 
negative. The area —d < x < 0 is in the normal state as 
before, A = C/g = 0, but the interlayer voltage is finite, 
E = Eo > 0. This region is a tunnel barrier between the 
superconductor and the domain wall at the interface to 
the third region x > 0, where A = Us = 0 and E = — Eq. 
In this situation guided modes exist at |£| < min(Eo, A) 
because states above Vb can propagate in the normal re¬ 
gions and states above A can propagate in the supercon¬ 
ductor. Because of the tunnel barrier the guided modes 
are only weakly affected by the superconductor and we 
can apply standard quasidegenerate perturbation theory 
[30l|3T], for which the unperturbed Hamiltonian Hq is ob¬ 
tained from i?BdG by setting A — U = 0 everywhere and 
so the perturbation Hi, which adds the missing parts, 
is finite only at x < —d. As a result of the perturba¬ 
tion the electron and hole states of the channel acquire 
a finite overlap A„„/(py), where n,n' label the subgap 
bands. This allows for particle number non-conserving 
processes, i.e.. Cooper pair transport. To first order only 
the electron and hole states belonging to the same subgap 
band mix, An = A 22 , A 12 = 0. This agrees with the 
result one expects when introducing superconductivity 
phenomenologically by constructing the BdG equation 
directly from the guided modes with a uniform pairing 
AiiTx. The second order corrections, which take into ac¬ 
count the modification of the wavefunctions due to the 
superconductor, however, reveal that the situation is dif¬ 
ferent in the geometry we consider. The electron hole 
overlap differs in both band^ An 7 ^ A 22 , and band mix¬ 
ing is finite, A 12 7 ^ 0 (Fig. 0 [52]- This is confirmed by 
the full dispersion relation of i?BdG we obtain by match- 
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FIG. 2. Induced intraband superconductivity An (A22) 
in the ID channel at the respective Fermi point py = 
±2-4^^ OL and induced interband superconductivity A 12 
at the band crossing py = Q. For illustrative purposes we 
choose the bulk superconducting gap A = Vb and the dop¬ 
ing Us = 10Vb (A < Vb is equally feasible). The amplitudes 
decay exponentially with the separation d between the super¬ 
conductor and the channel because the F > 0 region acts 
as a tunnel barrier. Inset: In the normal state dispersion 
(dashed) two different-sized gaps open at the Fermi energy 
because An yf A22, shown for d = Addition¬ 

ally, Ai 2 opens a gap at the electron-hole crossing at py = 0. 
This point contributes significantly to Cooper pair transport 
because compared to the Fermi points the normal density of 
states is higher and because the energy is larger such that 
the bound states extend further into the bulk, increasing the 
coupling to the superconductor. 


ing the 4-component spinor and its derivatives at both 
interfaces numerically (Fig. inset): two gaps of differ¬ 
ent size open at zero energy (An and A 22 ) and two gaps 
open at zero momentum where electron and hole states 
from different subgap bands cross (A 12 = A^i). This 
means that there is Cooper pair transport at zero energy 
as well as at the finite energies ±Vb/-\/2- 

Cooper pair transport .—We use Fermi’s golden rule to 
calculate the Cooper pair current I = 2e^j-(FF^ — 
Wy^Pi, where Wf. = ^\{f±\HT\i) p(5(e/ - ei) is the 
transition rate from an initial state i with probabil¬ 
ity Pi at energy et to the final state f± with 2 more 
(less) electrons at energy £f. The tunnel Hamiltonian 
Hp comprises the particle number non-conserving terms 
of the second-quantized perturbative model with elec¬ 
tron operators and hole operators = 

sc^^^ _y{—k), where k = Py, 

Ht= Y. '^nn'{k)scl^^yk)cY,-s{-k)+CL.c. (4) 

Xv'^’ci'ks 

Because the superconductor interface has a finite width 
w, we restrict the pairing amplitude in real space A(x, x') 
to x,x' S [—w/2,w/2]. In momentum space (sup¬ 
pressing all indices) this amounts to A(fc)cfcC_fe —> 


Afefc/CfeCfc' with A(fc) from the microscopic calcula¬ 
tion and 


Akk' = A 




sm 


(l-k)'. 


sm 


[l + k' 


(i-m 


il + k')k 


where L is the total length of the system, which 
does not enter the final results, and we have ex¬ 
ploited that the integrand is peaked around A: « Z « 
—k'. With this, the rates for removing (adding) a 

Cooper pair, |i) —>■ |/) = Cx^J],{k)c"_^Ysik') I*)? become 

where at low temperatures the occupation 
probability (n^yjk))^ = 1 - (ny^k)), « - 

£^^{k)) with Sp, the voltage applied between the super¬ 
conductor and the channel. Rewriting the sum over mo¬ 
menta as energy integrals, the current becomes 


I = 


32e ttL^ ^ 

p6fi 

1 ria 

dh K 


h (27r)2 ^ j 
nn' 

1 Leo 

-Sp 

de 

de 


^ \Ann*{£^ S')! . 


(5) 


The combination of energy conservation and approximate 
momentum conservation implies that the pair tunneling 
probability |A„„/(e, —e)P has a single peak as a function 
of £ for each pair n,n\ Injection into the same subgap 
band, n = n', happens near £g" = 0, and into different 
subbands, n y n', near = —e^ = Fq/v^ (Fig. [^). 
Linearizing the dispersion ([^ around these points [31) . 
£ = £ 0 " +hv™ (fc —fcg” ), the tunnel amplitude becomes 
Ann’ = A„„/(A:o”')sin[(£ - e]^'^'){w/hv'^'^')]/L{e - 
and we obtain the conductance 


G « 4Gn V Tn 


Sn,(Sp-er) + Sy,(Sp + er) , ( 6 ) 


where Go = 2e^jh is the conductance quantum, dw{£) = 
sin^ £w/(fiz;o" ) /{nwe'^) becomes the delta func¬ 


tion for re —>• 00 , and T„„/ = 27rw|A(A:Q”')p/(Zi?;o"') is 
the effective tunneling strength. Note that the conduc¬ 
tance grows with the length of the interface. This is 
in contrast to conventional Cooper pair splitters, which 
suffer from an exponential suppression in the spatial size. 
The reason is that here Cooper pairs are split kinemat¬ 
ically only after having tunneled locally into the chan¬ 
nel, a process which can happen simultaneously along the 
whole interface. The conductance contains a central zero- 
bias peak and two characteristic side peaks [Fig. [^a)], 
which arise because of the special subgap band structure 
and which correspond to the injection points marked in 
Fig. [3b). The peak height is proportional to the induced 
superconducting pairings. A factor of 4 arises due to the 
spin and valley degeneracy and a factor of 2 due to pair 
transport. 

Local Andreev reflection and Cooper pair splitting .—In 
Eq. (jij) the singlet nature of the injected Cooper pairs is 
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FIG. 3. (a) Subgap conductance G and Cooper pair current 
7 of a to = 2?)h^/aJVo long interface between the supercon¬ 
ductor and the ID channel at a distance d = “iy/WaJVQ with 
A = Us = Vo. The peak structure reflects simultaneous en¬ 
ergy and approximate momentum conservation. The oscilla¬ 
tions are caused by the sharp boundary of the superconduc¬ 
tor region and vanish if an exponential cutoff is used instead 
(dashed). (b) Interpretation of Cooper pair splitting in terms 
of Andreev processes. Incoming holes (open circles) filled up 
to the bias are either transmitted (T) or locally Andreev- 
reflected (LAR). Ordinary reflection and crossed Andreev re¬ 
flection (CAR) are zero by the valley chirality. A LAR process 
creates an outgoing electron (filled circle) on the same side and 
no outgoing hole on the opposite side, which corresponds to 
an electron of opposite spin, momentum and energy (dashed 
arrow). These two electrons are spin entangled (text). 


manifest. It is well established that Cooper pair splitting 
is closely related to CAR [531131] ■ This applies if the 
dominant process is ordinary reflection. In our device, 
the ID channel with a proximity-induced superconduct¬ 
ing region is a NSN junction, in which only transmission 
through the S-region with amplitude t(s) and local An¬ 
dreev reflection (an incoming quasiparticle in valley K± 
is rehected as an outgoing antiparticle with opposite ve¬ 
locity in valley Kz^) with amplitude r(e) are possible. 
We consider the general scattering problem with Hnite 
CAR and normal reflection in the Supplemental Mate¬ 
rial [3l|. When a voltage bias Sfi is applied between the 
superconductor and the channel to extract Cooper pairs, 
the incoming modes are hlled with holes up to + <5/r. 
Without the superconductor, all are transmitted and hll 
the outgoing modes up to Ep + Sfi, which is equivalent 
to a Fermi sea for electrons with the Fermi energy 

Ep — Sfj, [ 331133 ] 

n n c:(-e)i). (7) 

0<£<Sfj, 0<£<6fj, 

OLS OlS 

Here, |) is the quasiparticle vacuum with respect to the 
Fermi level Ep of the superconductor and /is^^^(e) = 


sc^{^{—s) creates outgoing holes with spin s at energy 
Ep -|- £ in the left/right lead, which is the same as an¬ 
nihilating an electron with opposite spin — s at energy 
Ep — £. We drop the valley index which is hxed by the 
requirement that outgoing modes move away from the 
superconducting region and the band index for simplic¬ 
ity [3T]. Due to the proximity effect LAR becomes hnite. 
The key observation is that when LAR occurs, no hole 
with spin s at energy Ep + e is transmitted to the other 
side. The outgoing mode is therefore occupied by a spin 
—s electron at energy Ep — e [Fig. |^b)]. To see this, we 
use Eq. Q to write the outgoing state in terms of 

m, 


n (sic“,(-£) -brc“1'(£))c“t(-e) 


0<£<5/J, 

OiS 

n 

0<£<SfJ, 


-I-r^C?’'’(£)c“’*’(-£)c“’*’(£)c?’'’(-£) 


+ 


r7(c“^(£)c“^(-£)-c“^(£)c“’^(-£)) l)^^. (8) 


If r is small, it becomes [1 -b ~ 

c;\e)cf{-e))+0{r^)]\)s^, where the desired nonlocal 
singlet state is explicit. This corresponds to a situation, 
where individual splitting events are well separated and 
it is meaningful to talk about pairs. In this regime of 
interest the perturbative result from the previous section 
holds. Only the emitted pairs contribute to the shot noise 
of the scattering state. In the opposite limit of perfect 
LAR with 0{t) ^ 0, 0(r) ^ 1, the outgoing state is a 
nonentangled product state. LAR is most pronounced 
at energies £ = 0 and e = iVo/v^ [Fig. |^a)] where 
the superconductor opens gaps A„„/ in the spectrum for 
the case of an inhnitely long (w —> oo) tunnel-junction 
(Fig. I). The LAR process becomes weak for all energies, 
when w falls below the coherence lengths Iiuq ” /A„„'. 

Conclusion .—Our setup allows for highly efficient cre¬ 
ation of nonlocal spin-entangled electrons without the 
need for repulsive interaction or energy hlters. We note 
that the topological channel can be created electrically 
in the bulk of the BG sample, completely avoiding sharp 
sample edges, the main source of intervalley scattering 
[36] . which could reduce the splitting efficiency. More¬ 
over, using an electrically tunable channel geometry bal¬ 
listic beamsplitters could be created to prove the spin 
entanglement via noise m, so far an elusive goal. The 
spin relaxation and decoherence in BG are expected to 
be weak due to the small spin-orbit coupling [35H3n] and 
the sparsity of nuclear spins. 

We thank A. Baumgartner, P. Samuelsson, 
C. Schonenberger, and A. Levy Yeyati for helpful 
discussions and acknowledge support from the EU-FP7 
Project SE2ND, No. 271554, the DFG, Grant No. RE 
2978/1-1 and Research Training Group GrK1952/l 
“Metrology for Gomplex Nanosystems”, and the Braun- 






















5 


schweig International Graduate School of Metrology 
B-IGSM. 


Supplemental Material 
Perturbation theory for the proximity effect 


Here we give the details of the calculation of the prox¬ 
imity induced amplitudes A. Because the superconduc¬ 
tor is only weakly coupled through the tunnel barrier 
—d < X < 0, we can derive an effective ID model via 
low-order quasi-degenerate perturbation theory. We split 
^BdG “ ^0 + Hi into two parts, where Hq is diagonal in 
the eigenbasis < 1 )°’"^ (x), and Hi = 0 (—x — d)[{—Us — 

Voaz)Tz + Axa;]. Since Hi is diagonal in spin and valley, 
we suppress the indices , s in the following. To first 
order in Hi, 

Hl2,{Py) = I dx^°f{x)Hiix)^°f{x) (Al) 
and to second order. 


H^3iPv) = 1^2 f 

X J dx^ll^^ix)Hiix)<i>lf (x), (A 2 ) 


where <I)p^_p^(x) are the unperturbed free states above 
the gap (|e| > Vq) with real and Py at energy = 

ija^{p^ + We impose the quantization condi¬ 

tion px = 27rn/L and normalize the extended wavefunc- 
tions according to 

quantization length L and the highest momentum px are 
increased until the second order matrix elements con¬ 
verge. To study Cooper pair transport only the parts of 
iy( 2 ) 

are relevant which are proportional to Tx, i.e., 
they mix electron and hole states and therefore change 


the particle number. The relevant momenta Py are close 
to the crossing of the respective electron and hole band 
(see the discussion on approximate momentum conserva¬ 
tion in the main text). This can involve one band. An = 

h\21{Pf) + H^21{pf) and A 22 = (-Pf)+ ^^24 (“Pf), 

or both, Ai 2 = A 21 = ( 0 ) -|- where ±pf are 

the Fermi points of the unperturbed dispersion, Eq. (§. 


Linearized subgap dispersion 

The linearized subgap dispersion, Eq. ([^, reads 

(Bl) 

around the Fermi points and 

F;«±^-x.2-i/4yi^Pp (B2) 

around zero momentum. The coefficients Eg" and Xg" 
used in the transport calculation, e.g., Eq. , can be 
read off immediately. 


Local and crossed Andreev reflection 

In the most general case the incoming holes in a NSN 
junction can be transmitted {thh), reflected {ruh), or un¬ 
dergo local {reh) or crossed (teh) Andreev reflection. The 
outgoing state is 

(thhKHe) + tehcfis) + rhhK^e) + re;,c“’l'(E)^ I) . 

escx. 

(Cl) 

Rewriting the hole operators hi in terms of electron op¬ 
erators c, and the Eermi sea |) in terms of the lowered 
Eermi sea |)^p as explained in Eq. Q in the main text, 
we arrive at 


n [i^hh+^lh) + (rlh - t2^)cfl(E)c((l(-E)cf'l(E)c(:t(-E) 

SS 

+s{rehrhh - thhteh){c’^\£)cl\{-e) + cf\£)c^\{-e)) 


+ s{tehrhh - thhreh)[cf\£)c^,.[-£) + cf'* (£)c^l(-£)) 



(G2) 


The first line contains the product state contributions, ~ 1, realized in Y-junction Cooper pair splitters, the 

the second line local pairs, and the third line nonlo¬ 
cal pairs. In the conventional reflection-dominated case, 
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leading order contributions are 

^ H {'^ehsic^\e)c^l{-s) + cf'f{s)c^l{-e)) 

SS 

+tehsicf\e)c^l{-e) + c^\e)c^l{-e)) 


/Sp. > 

(C3) 


i.e., LAR produces local pairs and CAR produces non¬ 
local pairs. In the transmission-dominated situation, 
thh ~ 1, the situation is reversed: the leading order is 


^ H {^ehs{c^He)c^l{-e) + cf1'(e)cf!^(-e)) 

SS 

+rehs{cf\e)c^l{-e) + c^\e)c^l{-e))^ 


I) 


6fi ’ 


(C4) 


so LAR produces nonlocal pairs and CAR produces lo¬ 
cal pairs. In the situation discussed in the main text, 
both CAR and reflection are forbidden, ruling out lo¬ 
cal pairs to all orders, as long as the valley symmetry is 
obeyed. Generally speaking it is undesirable to have si¬ 
multaneously strong ordinary reflection and strong LAR 
or to have simultaneously strong transmission and strong 


CAR to build a Cooper pair splitter useful to create spin 
entanglement. 


Multiband Andreev reflection 

The notation becomes more cumbersome, when both 
snbgap bands are considered but the considerations are 
completely analogous. Without superconductivity the 
outgoing scattering state is 

n Kiie)\)^ n c(£)i)^i)^, (Di) 

0<£<Sfi 0<£<6fi 

ocsn asn 

where n € {15 2} is the band index. In the presence of 
the superconductor, the transmitted holes can change the 
subgap band from n to m with an amplitude tmn(e)- Like 
in the one-band case, whenever the energy of an incoming 
electron is such that the spectrum of the S region has a 
gap, the transmission amplitude is exponentially 

suppressed with the length of the proximity region, and 
due to unitarity there is a finite amplitude r„m(£) for 
the spin-s hole to be Andreev reflected locally as a spin-s 
electron at energy Ep + s: 


I out) 


^ ' fs tmn C^,-s( E rmn \) 

0<e<Sp m 
arts 

0<€<Sfi m m' 

arts 

n [(^12^21 - ^11^22) + (riir22 - ri2r2i)c“J’(e)c2j'(e)c“l,(-£)c2,l,(-£) 

0<£<Sfl 

(DCS 




I) 


' 


(D2) 


Writing out the spin part of the product explicitly, the 
nonlocal singlet nature of the injected Cooper pairs be¬ 
comes obvious: 


|0Ut) — (^12^21 — ^ 11 ^ 22 )^ 

— 5^<e<5ii 

+ (tl2t21 - tiit22)^{-lT{rnltrn2 - (e)c^\ (-e) - 


0{r^) I),,. 


(D3) 


The higher order terms in r contain multiple Cooper pairs and are not necessarily entangled, e.g., the 0{r^) contri- 
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bution is a pure product state in which all states in the 
left/right lead at energy Ep ± e are occupied. 
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